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Abstract 

We extend the Wolff algorithm to include correlated spin interactions 
in magnetic systems. This algorithm is applied to study the site-bond- 
correlated Ising model on a two dimensional square lattice. We use a finite 
size scaling procedure to obtain the phase diagram in the temperature- 
concentration (T c , C) space. 
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For randomly diluted magnetic systems the critical concentration where the 
magnetic order vanishes is of great interest. The dependence of the critical 
temperature with the concentration is ruled by the topological properties of the 
lattice and by the symmetry of the interaction Hamiltonian. 

In the simplest description of a dilute magnet, the occupied lattice sites 
correspond to the magnetic atoms and the empty sites are associated with the 
presence of nonmagnetic atoms (impurities). Interactions are short ranged usu- 
ally of the exchange interaction type. The strength of the exchange interaction is 
not affected by their neighbors. At T = 0, the cluster structure is well described 
by the ordinary site percolation model. 

However, it has been found experimentally that sometimes the local environ- 
ment may modify the exchange coupling constant between two atoms or may 
even supress their magnetic manifestations. This means that interations involv- 
ing more than two atoms are present in the system. To take these effects into 
account, two kinds of correlations were proposed: the long and the short ranged 
correlation models. The bootstrap percolation model Q is a good example for 
the former and the site-bond-correlated model j^, || for the latter. Long-ranged 
correlation may change the critical exponents or even the order of the transition 
of its corresponding uncorrelated model. In contrast, the short-ranged corre- 
lation seems to be unable to such drastic effects. Anyway, for these systems, 
dilution as well as correlation play an important role. 

An analysis of 19 F NMR linewidths in the randomly diluted magnetic 
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system KNi x Mg\- x F^ and in the isostructural compound KMn x M g\- x F^ ( 
where x is the concentration ) show remarkable differences in their proper- 
ties |j. For instance, the concentration or percolation threshold where the 
magnetic order ceases are different. Further, the former displays an upward cur- 
vature in the temperature-concentration plane which is absent in the latter. In 
order to explain these facts, Aguiar et al. [^) proposed a dilution model where 
the exchange coupling constant of two ions Ni 2+ depends directionally on the 
magnetic attributes of theirs nearest neighbor atoms. That is the magnetic spins 
are correlated. In a subsequent work |^], a parameter a was introduced which is 
a measure of the correlation strength. When a > or a < we have ferro or an- 
tiferromagnetic correlation, respectively. Due to the presence of both ferro and 
antiferromagnetic coupling constants, the antiferromagnetic correlation brings 
competition to the system and a spin glass behavior is expected. 

The thermal properties of the site-bond-correlated Ising model (hereafter de- 
noted SBC) was investigated using several approaches: mean field ||, Honmura- 
Kaneyoshi effective- field || , renormalization group |(| Q and Monte Carlo renor- 
malization group ||] . Any of these techniques lead to almost the same physical 
scenario with the exception of the SBC defined on the Bethe lattice ||. Here 
there are some unexpected features which seem to be due to the pathological 
geometry of the Bethe lattice. Also the SBC cluster characteristics has been 
studied and its connection with the percolation problem. For the square lattice, 
it is now well established that there are two kinds of percolations: the usual site 
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percolation with threshold p c ~ 0.592 when < a < 1 and a correlated percola- 
tion with threshold p c ~ 0.740 at a = |l(|. There are now strong indications 
that both percolations belong to the same universality class [[0], [l2| . 

As far as we know, reference [ fl3| is the only work which treats the thermal 
properties of the SBC model by using the Monte Carlo technique. In that paper, 
the parameter a is resctricted to be null. Everyone knows that is very difficult 
to simulate near phase transitions due to the emergence of the critical slowing 
down phenomenon. In the SBC this becomes even worst since dilution together 
with correlation conspire to weaken the Metropolis technique. Fortunately, for 
the so called cluster algorithms the dilution turns these algorithms even more 
robust |l4j]. In the Ising model (diluted or not) the Metropolis algorithm is 
a local Monte Carlo method where only one spin is flipped each time. For 
the cluster algorithms however, the entire cluster can be flipped. Indeed, there 
are actually two kinds of cluster algorithms: the single-cluster algorithm like 
the Wolff [jTB) algorithm and those involving multiple clusters as the Swendsen- 
Wang ]l6| and the invaded cluster jl7j algorithms. In this paper we generalize 
the Wolff algorithm to include correlation. We use our results to construct the 
SBC phase diagram for many concentrations and for the correlation a in the 
interval [0, 1]. 

In spite of the great interest received during the last decades by the critical 
slowing down phenomenon, it is still a hard question to find under what con- 
ditions it can be reduced. Two independent strand of thought had contribute 
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to enlighten the problem. The first was the search for the so called "physical 
cluster" of the Ising model. The idea is to find a percolation problem in abso- 
lute harmony with the thermal problem: the temperature where the percolation 
threshold occurs (T p ) must coincide with the critical temperature (T c ) and also 
their corresponding critical exponents need to be the same. The "geometrical 
cluster", i. e., those clusters formed by grouping all nearest neighbors up (or 
down) spins were discarded since in 3-d they neither have T p = T c nor are their 
exponents the same. The first explicit construction of "physical clusters" was 
proposed by Coniglio and Klein jl8| . Parallel to this research line, Fortuin and 
Kasteleyn fill introduced the random cluster model. They proved that the sus- 
ceptibility of the Ising model is equal to the mean cluster size of the random 
cluster model. Gathering these ideas, Swendsen and Wang (T(| developed a new 
algorithm for the Potts model. In this algorithm, clusters of spins in the same 
states are grown and flipped. Thus, a number of spins are update in a single 
move and the correlation time is reduced. 

In 1989, a single- cluster Monte Carlo algorithm was introduced by Wolff 
i jTs) for the O(N) spin models as a variation on the Swendsen- Wang scheme 
for the Ising model and has proved to be even more effective. In the Wolff 
algorithm only one cluster is formed and flipped with probability 1, whereas in 
the Swendsen- Wang dynamics all percolation clusters are formed and flipped 
with probability 1/2. According to Tamayo et al. [^0] the main reason for a 
better performance of the Wolff algorithm is that the mean size of the clusters 
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flipped is significantly larger than in Swendscn-Wang case. 

In the SBC model the presence of nonmagnetic impurities in the neighbor- 
hood of a given pair of nearest-neighbor magnetic atoms can modify the strength 
of the exchange interaction between the two atoms. Moreover, in the limit of 
strong correlation a, the correlation can even suppress the relevant exchange 
interaction. 

The model Hamiltonian is the following 



i,s 

where Oi = ±1 and 6 denotes an elementary lattice vector. The exchange 
interaction Jis + s is given by 



where J > 0. The random variables £j can take values one with probability 
C and zero with probability 1 — C, where C is the concentration of magnetic 
atoms. The parameter a correlates the interaction between sites i and i + S 
with the magnetic occupancy of the sites i — 5 and i + 25. The uncorrelated 
dilute Ising model is re-obtained in the limit a = 1. For < a < 1, the bond 
between i and i + S is only weakened by the absence of a magnetic atom at 
i — S or i + 25. The limit a = corresponds to the maximum correlation, i. e., 
two magnetic first neighbor sites are connected by an active bond only if their 
nearest-neighbor sites along the line joining them arc also present. 




(1) 



Ji,i+6 = Js%£i+S [(1 - a) £i-S£i+2S + ct] 



(2) 
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In the Wolff algorithm for the pure Ising model, first a site is randomly 
chosen in the lattice. A nearest neighbour site will be added to the cluster with 
an activation probability p = 1 — e~ 2K (where K = k J bT ) if it is in the same 
spin state. This procedure is repeated until no more sites can be incorporated 
to the cluster. The whole cluster is then flipped. Following Fortuin-Kasteleyn 
[ fl9| , we derive the bond activation probability: 

!Pi.i+s = if <Ji ^ a i+ 8 

Pi.i+s = 1 - e~ 2K if ^ = <r i+s and £i-se i+2 s = 1 
Pi.i+s = 1 - e~ 2aK if <7j = a l+s and e^s£i+2S = 

between sites i and i + S. The definition of these activation probabilities guar- 
antees correctly the applicability of the cluster dynamics and constitutes the 
Fortuin-Kasteleyn mapping for the SBC model. 

We simulate the SBC model for various values of correlation (a) and con- 
centration (C) on a square lattice of size L — 50, 100, 150 and 300. A ran- 
dom uniformly distributed magnetic sites configuration (or simply sample) is 
generated with an occupation probability C. Over this quenched (geometric) 
configuration, an initial spin configuration is chosen with half of the spins up. 
One Wolff's cluster is then constructed using the activation probabilities just 
described. The entire cluster is then flipped and the new magnetization is deter- 
mined. We call the sequence: cluster construction + flipping + magnetization 
measure one iteration or Monte Carlo step. Our first 1000 iterations were dis- 
carded waiting the system to achieve the thermal equilibrium. The remaining 
iterations were used to measure the mean magnetization and its fluctuation - 
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the magnetic susceptibility. On these quantities, another average over the sam- 
ples was necessary in order to anneal the geometric influence (always present 
in diluted systems). Finally, we scan the coupling K to find K max where the 
susceptibility is maximum. In Table 1 we show the number of iterations and 
realizations necessaries to reduce the error of K max to ~ 0.01. Of course, they 
depend on the values of a and C. 

Table 1 to be inserted 

The determination of K max was done for fixed values of a, C and L and then 
extrapolated to the thermodynamic limit L — > oo through the BST algorithm 
[ pT| . The BST is a useful algorithm to extrapolate physical quantities that 
converge obeying a power law F(L) — F(L = oo) + AL~®. It allows a reliable 
determination of critical parameters in the limit L — > oo and its versatility 
becomes more pronounced if there are only very short sequences available. 

Figure 1 to be inserted 

The phase diagram is shown in Fig. 1. The initial slope (l/T c )(dT c /dC) at 
the pure C = 1 case increases with the degree of correlation, which is supported 
by both experimental M and theoretical || ||, |l3) data. At T = 0, we obtain 
two distinct percolation threshold: one for a = and another for ci / 0. It 
can be seen in the figure that the curves for intermediate correlation < a < 1 
have an upward curvature ( not present in extreme the cases a — and a = 1) 
in agreement with the experimental results for the KNi x Mgi- x F^ compound 
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0. In summary, we have generalized the Wolff's algorithm for systems where 
both dilution and correlation are present. This allowed us to get the phase 
diagram of the site-bond-correlated Ising model on a two dimensional square 
lattice. The natural question on how the new algorithm changes the dynamical 
critical exponent and the autocorrelation time is now under investigation. 

We acknowledge Conselho Nacional de Desenvolvimento Cienti'fico e Tec- 
nologico (CNPq) and Fundagao de Amparo a Pesquisa do Estado de Sao Paulo 
(FAPESP) for financial support. 
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Table Caption 



Table 1 - Range of the number of samples N s and number of iterations 
Ni used in this work. The numbers N s and Ni augment for smalls a (strong 
correlation) and C (strong dilution). 
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Figure Caption 



Fig. 1 - The SBC Ising model phase diagram for many values of a. The 
percolations threshold are 0.740 and 0.592 for a — and < a < 1, respectively. 
The dashed lines are only a guide to the eye. 
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L 


N s 


N< 


50 


20 - 150 


3000 - 5000 


100 


05 - 50 


3000 - 5000 


150 


05 - 50 


3000 - 5000 


300 


05 - 50 


3000 - 10000 



Table 1 
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Figure 1 
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Figure 2 
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Figure 3 
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